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ABSTRACT 

We  study  the  class  of  functions  < x < “>)  having  the  following 

properties: 

(i)  The  restriction  of  S fx)  to  every  unit  interval  (V,v+1) (ve  Z)  is  a 

n 

polynomial  S (x)  of  degree  not  exceeding  n. 
n,  V 

(ii)  If  h > 0 and  n h < 1,  and  S ,(x),  S , , (x)  are  two  successive 

ri/v  n^v+i 

polynomials,  then  they  satisfy  the  Abel  interpolatory  conditions 

s'  ,(v  + 1 + kh)  = s'  ' (V  + 1 + kh)  (k  = 0,l,...,n-l). 
n,v  n,v+l 

These  are  a kind  of  discrete  splines  that  may  be  called  cardinal  Abel  splines 
It  is  shown  that  the  theory  of  the  usual  cardinal  polynomial  splines  of  degree  n 
(See  [6])  carries  over  to  the  new  class  of  cardinal  Abel  Splines.  The  latter  clas 
reduces  to  the  former  if  we  let  h 0.  We  obtain  in  this  way  generalizations  of 
the  exponential  Euler  polynomials  of  the  exponential  Euler  splines,  eigensplines 
a. s. f . 
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DISCRETE  EXPONENTIAL  ABEL-EULER  SPLINES 


R.N.  Mohapatra  and  A.  Sharma 

1.  Introduction.  Discrete  splines  made  their  appearance  first  in  the 
literature  when  Mangasarian  and  Schumaker  [2]  brought  out  their  minimization 
properties  and  their  relations  with  summation  formulae.  Later  Tom  Lyche  [1] 
made  a systematic  study  of  the  subject  in  his  dissertation.  More  recently 
Tzirabalario  [8]  has  studied  discrete  cardinal  splines  and  has  brought 
out  an  analogue  of  Schoenberg's  exponential  Euler  splines.  He  has  shown  that 
the  methods  of  Schoenberg  can  be  adapted  to  solve  the  interpolation  problem 
to  data  of  power  growth  by  discrete  splines  with  the  same  po^^^er  growth  and 
having  integer  knots. 

One  of  the  significant  features  of  the  discrete  splines  of  degree  n so  far 
available  is  that  two  polynomial  components  of  a discrete  spline  of  degree  n on  two 
adjacent  intervals  determined  by  three  successive  knots,  say  x^, 

agree  in  n equidistant  points  x^  + kh,  k = 0,l,***,n-l.  This  require- 
ment is  equivalent  to  requiring  the  first  n differences  of  the  polynomials 

on  these  n points  to  vanish.  Thus  the  two  polynomial  components  differ 

n-1 

only  by  an  expression  c H (x  - x - kh)  where  c is  a constant. 

k=o  '' 

It  appears  that  there  is  no  special  merit  attached  to  this  particular 
constraint.  The  object  of  this  paper  is  to  replace  this  requirement  by  another 
one  related  to  the  Abel  interpolation  problem.  It  turns  out  that  the  analogues 
of  Schoenberg's  result  on  Cardinal  splines  come  out  equally  easily  in  this  case 
as  well. 


Ku  fha  IlnlfcAd  StAtAfl. 
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It  also  follows  from  this  treatment  that  there  can  be  many  discrete 
splines  depending  upon  the  nature  of  Interpolatory  restraints  which  weld 
the  two  adjacent  polynomial  components.  The  discrete  splines  considered 
by  Mangasarian  and  Schumaker  [2]  and  Lyche  [1]  give  one  example  of  this 
observation. 

In  Section  2,  we  give  the  preliminaries  and  formulate  two  problems  of 

interpolation.  We  also  define  the  B-splines  Q^,,(x)  which  form  a basis  for 

n+i 

the  discrete  splines  studied  here.  Section  3 deals  with  exponential  Abel- 
Euler  polynomials  and  their  properties.  We  introduce  the  polynomials 
^n  h/n^*’^^  show  that  for  h < x < 1,  these  polynomials  have  real 

simple  negative  zeros.  In  Section  4 we  use  the  results  of  Section  3 to 
define  the  exponential  Abel-Euler  spline  and  study  its  convergence  as  its 
degree  tends  to  infinity.  In  Section  5 we  sketch  the  outlines  of  a solution 
of  Problem  A of  Section  2.  We  define  t-perfect  discrete  Abel  splines  in 
Section  6 in  analogy  with  the  study  of  Sharma  and  Tzimbalario  [7].  We  propose 
the  general  problem  and  solve  it  for  r = n - 1.  For  r ^ n - 2,  we  are 
not  yet  able  to  resolve  the  problem. 

2.  Preliminaries.  Let  denote  the  set  of  polynomials  of  degree  £ n 

and  let  P denote  the  class  of  functions  S (x)  whose  restriction  S (x) 
n n n,  V 

to  (v,v+l),  for  all  integers  v,  belongs  to  tt  . For  a given 

n 

real  positive  number  h with  nh  < 1,  we  denote  by  the  class  of  functions 

n 

S^(x)  which  satisfy  the  following  conditions: 

(1)  S (x)  c P , 
n n 

and  for  v = 0,  ±1,  ±2,***, 


I 
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(ii)  (v+l+kh)  = (v+l+kh),  k = 0,l,-“,n-l, 


where  S (x)  and  S ,,(x)  denote  the  restrictions  of  S (x)  to  (v,v+l) 
n,v  n,v+l  n ’ 


and  (v+l,v+2)  respectively. 


We  shall  call  S the  class  of  Cardinal  Abel-discrete  splines.  As 
n 


h ->  0,  reduces  to  the  class  of  Cardinal  splines  in  the  sense  of  Schoenberg 


Following  Schoenberg  [6],  we  propose 


Problem  A.  Let  0£a<l  and  let  y = {y^},  v = 0,±1,«**,  he  an 


arbitrary  infinite  sequence 


(2.1) 


|yj  = Odvl^). 


Find  a function  S (x)  e S such  that 
n n 


(2.2) 


S (v+a)  = y , V = 0,±1,*«*, 
n V 


(2.3) 


ls^(x)|  = 0(|xd)  as  X -»■  ±«>, 


If  Yy  = ^ » A ^ 0,1,  then  the  data  Is  not  of  power  growth  and  leads 


to  Problem  B. 


Problem  B.  For.  a given  o,  0 < a < 1,  find  a function  S (x)  e S"  which 

n n 


satisfies  the  conditions: 


(2.4) 


S^(v+a)  = X , V = 0,±1,' 


(2.5) 


S (x+1)  = X S (x). 
n n 


We  shall  first  turn  to  the  solution  of  Problem  B which  in  turn  will  lead 


to  the  solution  of  Problem  A.  We  shall  require  the  fundamental  polynomials  of 


Abel-Gontcharoff  interpolation: 


“ x(x-nh)"~^  (n  = 2,3,***). 

It  is  easy  to  verify  that 
(k) 

(2.6a)  kin. 

Set 


(2.7) 

where 


J*0  »»>»•  ^ 


(X  < 0), 


It  is  easy  to  see  that  Qj^^j(x)  c s|j  and  has  the  support  [0,n+l].  Following 
Schoenberg  (6]  it  can  be  shown  that  every  S(x)  e has  a unique  represen- 
tation in  terms  of  qJJ^j^(x)  and  its  translates,  l.e., 

(2.8)  S(x)  - X c Q^2<x-j). 

In  view  of  (2.8),  the  Problem  B can  be  easily  solved.  Indeed 


(2.9)  S_(x)  - c 


<•  jL 


satisfies  the  relation  (2.5),  and  (2.4)  will  be  satisfied  If 

(2.10)  c^  I qJJ+i(o-J)  - 1. 


We  shall  show  in  Section  4 that  for  nh  < a < 1, 


1 


-5- 


1 

I 


P 


, _ I 

t I 

k 


I 


i 


However  the  representation  (2.9)  is  not  satisfactory  for  the  study  of  the 
convergence  problem.  We  shall  obtain  another  representation  in  the  next  section. 


3.  The  Exponential  Abel-Euler  Polynomials.  We  shall  consider  the  monic 
polynomials  ^(x,X)  of  degree  n in  x which  are  determined  by  the  con- 
ditions: 

(3.1)  A^*^i(l  + kh,  X ) = X A^*^i(kh,  X ) (k  = 0,l.*-*,n-l). 

n,n  n,n 

We  shall  call  these  polynomials  exponential  Abel-Euler  polynomials.  Since  the 

polynomial  G , (x)  is  of  exact  degree  n,  the  polynomial  A . (x.X)  is 
n>n  n,h 

represented  by 


n,n  n,n  'I'  1 n-l,n  n-i  l,h  n o,h 


The  conditions  (3.1)  lead  to  the  following  system  of  equations  in  the  c^'s  on 


using  (2.6a): 


(3.2)  G^^h^^^'^^3‘'lVl,h^^^'^‘"'^Cl^Vl^l,h^^^  + % “ “ 1,2." -.n), 


It  follows  from  (3.2)  that  if  X * 1,  the  c 's  are  uniquely  determined. 


In  order  to  obtain  some  simple  properties  of  the  polynomials  A , (x.X)  we 

n.h 


observe  that  these  polynomials  are  given  by  the  generating  function: 
(3.3) 


k -1  "“A  . (x.X)  . 

X - 1 xz  r n.h  _n  _nzh 


X - e 


= I 

n=o 


n! 


z e 


This  follows  from  (3.2)  on  observing  that  the  sequence  {c  } is  independent,  of 


n and  hence  (3.2)  gives 

‘ I ^ ^ ■ I ^ ^ j.  Q 


I ; 
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Setting  c = 1, 
o ’ 


f(t)  = j;  c , 

V=o  ''' 


«<'>  = I <=„>«>  ^ • 


we  get 


(3.4)  f(t)  = 


A - 1 
A-g(t) 


zh 


Putting  t - ze  in  (3.4)  we  easily  get 


(3.4a) 


n nzh 


X-1 


X-e"  n=o  " 


= I - " " 


It  is  known  and  is  also  easy  to  verify  that 


(3.5) 


e = 1 + xze  + 


...  . x(x-nh)"'^  n nzh 
• • • + : z e + • • • 


n 


(|zh|  < 1). 


Formula  (3.3)  now  follows  on  multiplying  (3.4a)  and  (3.5)  and  using  Cauchy 
product  on  the  right 

We  formulate  some  properties  of  A . (x,X)  in 


THEOREM  1.  The  polynomials  A^^^(x.X)  satisfy  the  following  relations 


(3.6) 

An  „(K.X)  . 

(1-X)  J]  a'*  G (x-v-1), 
v=o 

(3.7) 

*11,11^*’  ' 

= (-1)"  A .(l-x,X), 
n,-n 

(3.8) 

Ani^d+x.X) 

- X An  ,,(x.X)  - (l-X)Gn  ,,(x). 

(3.9) 

”x 

■ " Vl.h'*-*'- 

p*aMMfn 
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Proof:  We  observe  that  from  (3.3)  we  have  on  using  (3.5), 


n=o  * 


nzh  , , ^ . 

^ 6 _ 1 ^ g(x-l)z 


1 - Ae 


-z 


= (1-X)  I A^ 
v=o 


= (1-A) 


I 

v=o 


n=o 


n nzh 
z e 

n! 


= (1-A)  j; 

n=o 


n nzh 

z e r . V 

n!  Jo 


This  proves  (3.6). 


(3.7)  is  immediate  from  (3.3). 
From  the  identity 


A - 1 . g(x+l)z 

1 z 
A - e 


-(A-De’'"'  + M,\-l,)  gxz 
X “ 


we  get  (3.8)  on  using  (3.3)  and  (3.5). 

^ Lastly  (3.9)  also  immediately  follows  from  (3.3)  by  easy  manipulation. 
We  now  set 


(3.10)  “ (A-1)"  j^(x,A). 

U follovs  fro»  (3.3)  th^t  Is  a polynomial  of  degree  n in  1 

that 


(-1)"  G (x-1). 
n,n 


and 


In  particular  we  have 


Hq  ^ = Ax  + (1-x) 

, (x,A)  = A^  x(x-2h)  - A{2x(x-2h)  - 2x  - 1 + 2h} 


+ (x-1) (x-l-2h) . 


If  x = 0,  polynomial  of  degree  n-1  In  A. 

From  Theorem  1 we  can  now  get 


COROLLARY  1.  The  following  representations  for  n , (x,A)  ore  valid: 

n,h 

CD 

(3.12)  n (x.A)  = (-l)^(l-A)"'*'^  I G ^(x-v-1). 
n j n n * n 

v=o  * 


(3.13)  n (x.A)  = n!  J A^"  q"!*,  (u+l-x) . 
n , h n+1 

V=o 


Relations  (3.12)  follow  immediately  from  (3.6).  It  follows  from  (2.7) 


Hence  from  (3.12),  we  have 


1=0  ^ v=o 


I I (0„ 

1=0  ^ ^ 


which  easily  reduces  to  (3.13). 


We  now  formulate 


THEOREM  2.  The  polynomial  II  j^(x,A)  for  a fixed  x,  h < x < 1,  as 

" * n 

polynomial  in  A has  n real  negative  simple  zeros. 


Proof ; The  proof  of  this  result  depends  upon  the  following  relation: 


(3.14)  ^ n (x,X)  = X(1-X)D  II  . (x.X)  + {(n+l)X  - (1-X) (x-l-nh) }n 

n+1,  ~ ^ 

To  prove  (3.14),  we  differentiate  (3.12)  with  respect  to  X so  that 
<3.Ua)  . 


where 


,vn+l  ” 

l"  I vX''g  Jx-v-1). 

A n • n 

v=o  ’ 


Writing  vX  | (x-'l-nh)  — (x-v-l-'nh)  }x  in  and  observing  that 


we  have 


(x-nh)G  (x)  = G ( 

n , h , , nh 


X (x-l-nh)  n^^^(x,X)  - n 


From  this  (3.14a)  we  get  (3.14), 


Replacing  h by  — , we  have  from  (3.14), 


X(1-X)D^  n j^(x,X)  + {(n+l)X  - (l-X)(x-l-h)}n  j^(x,X) 

n»  ~ n,  — 


(x,X) 


We  observe  that 


(3.15) 


sgn  n ^ (x,0)  = 1, 


^1 


(x,-A)  = (-1)"  when  x > h,  (A  > 0 large). 


-10- 


If  A <X  <***<A<0  denote  the  negative  real  simple  zeros  of 
n n-i  i ^ ir 


n j^(x,A),  then 


"•n 


From  this  the  result  follows  by  induction  on  n on  using  (3.15). 


Remark.  Corresponding  to  each  zero  A^  of  II  ^(a,A)  for  h < a < 1,  we 

^>n 

have  a discrete  Euler  spline  S^(x)  which  will  satisfy  the  conditions: 


S^(a+k)  =0,  k = 0,±1,±2,-*-, 

and 


S (x+1)  = A S (x). 

V V V 


These  discrete  Euler  splines  have  exponential  growth  at  +«>  or  which 

depends  upon  |a^|  > 1 or  <1.  These  n discrete  Euler  splines  form  a 

basis  of  the  class  S°  , of  discrete  splines  where 

n,h  ^ 

S°  = {s(x)|s(x)  e S^,  S(a+k)  =0,  k = 0,±1,***}. 

Ti  j n n 


4. 

spline 


Exponential  Abel-Euler  Spline.  We  consider  the  discrete  exponential 


I Cl'-''' 


\)=—C» 


which  satisfies  the  functional  equation 


(4.1) 


<l>  , (x+l,A)  = X <J>  , (x,A). 

n,n  n,h  ’ 
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The  restriction  of  4>  , (x,X)  to  [0,11  is  A . (x,X)  and  « ,(x,X)  c 

_ £L  n,n  n 

"’n 

By  Theorem  2,  if  h £ a < 1,  X ?=  1,  then  A , (a,X)  * 0 if  \ \ , 

n V * 

"•n 

^ =-i. 2, which  are  the  zeros  of  n j^(a,X).  This  also  proves  (2.11)  in  view 
of  (3.10)  and  (3.13). 

Set 


(4.2) 


S , (x.X) 
n,h 


which  belongs  to  the  class  and  Interpolates  x’'  at  x = a ± v,  v 

We  shall  prove 


0,1. ■ 


THEOREM  3.  If  X is  not  a negative  nwnhev,  h < a < l^  and  S (x,X) 

n , h 

is  given  by  (4.2),  then 

(4.3)  Urn  S Jx.X)  = )X|^~“  eiu(x-a) 

n «,  n.h 

where  X = )x|e^'^. 


For  the  proof  of  Theorem  3,  we  shall  need  a relatioi 

nomials  A^  j^(x,X)  and  the  exponential  Euler  polynomials 

[6]  and  the  Fourier  series  representation  of  X'’^  A (x,X) 

n 


(4.4) 


From  (3.3),  it  follows  that 
” A (x.X) 


X - 1 xz 
e 

X - e^ 


= I 

n=o 


= I 


“ A , (x.X) 

n , h n 


n=o 


n! 


between  the  poly- 
A^(x,X)  of  Schficnherg 
which  is  1-pcriodic. 

nzh 


Hence  it  is  easy  to  see  that 


(4.5) 


n=o 


Inverting  this  relation  (Riordan  [4])  we  have, 

(4-6)  ^(x,X)  = J (”)(-!)"  ™ h”  ™ m n”  A^(x,X), 


m=l 


'm-' 


so  that 


(4.7)  A j^(x,X) 


m-1 


It  is  known  (Schoenberg  ([5],  p.  399  (7.14))  that 

n+i  * 0 1 ^ 1 1. 


, ^ “ 2irkix 

(4.8)  A (x.X)  = n!  ^ I ^ 

n A 


x: 


where 


iu 


X = log  T + iu  + 2nik  (X  = xe  , -x  < u < n) 


and  \ * 1.  Then  for  0 £ x _<  1,  it  follows  from  (4.2),  (4.7)  and  (4.8),  by 
easy  manipulation  that 

X -jn+l 


n,h  “> 


® . (A  '\n+l 

1 + I ^ B, 

—00 


kn 


1+  I 


2irkia 


X_'\n+1 


(^) 


kn 


where 


n-1 


W-”  »K„  • P JhC)  • 


n - m m 
X 
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and  ^ denotes  summation  for  all  k 0.  It  can  be  shown  with  some  elementary 

but  tedious  manipulations  which  we  omit  that  are  uniformly  bounded. 

Since  jx  /X  | <1,  it  follows  that 
' o k 

li»  S Jx.X)  . 1x1*-“ 


where  X = |x|e  This  completes  the  proof  of  Theorem  3. 


5.  Solution  of  Problem  A.  The  solution  of  Problem  A is  now  given  by 

THEOREM  A.  Let  h ^ a < 1 such  that  11^  j^(o,-l)  *■  0.  Then  for  every  set 

of  data  {y  } of  power  growth  there  is  a unique  function  S(x)  t s''  such  that 
V n 

(5.1)  S(\H-a)  = y^,  v = 0,tl,±2,*«*, 

(5.2)  S(x)  = O(Ixl^),  X ->■  +“ 


where 


y^l  = 0(|vr)  as  V ->■  ±« 


for  some  appropriate  y ^ 0. 


The  proof  of  the  theorem  follows  the  same  lines  as  that  of  the  corres- 
ponding result  of  Schoenberg  ([6],  p.  34)  and  will  be  omitted.  It  depends 

on  the  fundamental  discrete  spline  L(x)  c such  that 

n 


(5.3)  L(v+a)  =6  » v = 0,±1,±2,»«*; 

(5.4)  |l(x)1  < Ae"®!’"',  x€». 


|I,P]IL  I.  U , J I . 
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for  some  positive  constants  A and  B. 

It  Is  easy  to  see,  following  the  ideas  of  Micchelli  [3]  that 

i(x-a)u 


L(x) 


= j-  r 

2tt 

-TT 


Q(x,u)du, 


where 


2nikx 


J2(x,u)  = 


k=-“>  (l(u+2kTi)) 


n+1 


P^_l  [-ih(u+2kTr)) 


2TTika 


k=-“  (l(u+2k7r)) 


and  P , (x)  is  given  by  (4.9). 
n— i 


6.  t-Perfect  Discrete  Abel  Splines.  Let  m be  a natural  number  and  let 
r = For  a given  real  number  t we  shall  denote  by 


m,h,t 


(S(x)}, 


the  class  of  all  functions  S(x)  with  the  following  properties: 

(6.1)  The  restriction  ^(x)  of  S(x)  to  (v,v+l)  represents  a poly- 
nomial of  degree  m with  the  highest  term  t'*x™, 

( v+l+kh)  , k = 0 , 1 . • • . , r-1 . 
where  h is  a positive  number  with  rh  < 1. 


We  now  propose  the  following  problem: 


5- 


Problem  C.  Determine  S(x)  c S"  , ^ havina  the  least  t-norm  ||sl!  where 

ni^n^L  ^ tj'* 


I S lit  CO  = sup  - 

X e m ' t 


S(x) 
Jx]  • 


The  analogous  problem  for  h = 0 was  proposed  and  solved  in  [7].  When 
h ^ 0,  the  easiest  case  which  can  be  solved  is  when  r = m.  In  this  case 
we  can  easily  prove  following  the  argument  of  [7]  mutatis  mutandis 

THEOREM  5.  Suppose  m ^ 1 and  t is  a real  number  * 1 and  h a given 
real  number  (mh  < 1) . Then  there  is  a unique  discrete  Abel  spline 

ik  ^ ^ 

F (x)  € S , ^ which  minimizes  the  t^norm  wher*e 
m^n^t 

"k 

F (x)  = A.  j^(x,t)  , 0 1 X 1, 

m 

F*(x+1)  = tF*(x), 

and  given  by  (3.3)  with  X replaced  by  t. 

For  0 _<  r £ m-1 , the  problem  is  more  difficult.  We  shall  return  to 


this  problem  later. 
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